Kondo screening cloud in a one dimensional wire: Numerical renormalization group 

study 



o 
o 

(N 
S3 

a 
in 

(N 



Laszlo Borda 

Department of Theoretical Physics and Research Group "Theory of Condensed Matter" of the Hungarian Academy of Sciences, 
Budapest University of Technology and Economics, Budafoki ut 8. H-1521 Budapest, Hungary 

(Dated: February 5, 2008) 

We study the Kondo model -a magnetic impurity coupled to a one dimensional wire via exchange 
coupling- by using Wilson's numerical renormalization group (NRG) technique. By applying an 
approach similar to which was used to compute the two impurity problem we managed to improve 
the bad spatial resolution of the numerical renormalization group method. In this way we have 
calculated the impurity spin - conduction electron spin correlation function which is a measure of 
the Kondo compensation cloud whose existence has been a long standing problem in solid state 
physics. We also present results on the temperature dependence of the Kondo correlations. 

PACS numbers: 72.15.Qm, 73.63.Kv 
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Introduction — As being one of the most interesting 
quantum impurity problems, the Kondo effect [lj -when 
a localized impurity spin interacts with the itinerant elec- 
trons via spin exchange coupling- has been studied for 
decades both theoretically and experimentally 2] and has 
recently come to its renaissance when it has been pre- 
dicted to appear and observed Q in quantum dot sys- 
tems. Even though the Kondo effect seems to be well- 
understood now, still, there is a controversial question 
which has not been clarified yet: Whether there exists 
an extended "Kondo screening cloud" or not? 

In case of the single channel Kondo problem -when the 
local spin is coupled to one band of conduction electrons- 
the interaction with the electrons lifts the ground state 
degeneracy of the isolated spin resulting in an effective 
"screening" of the impurity spin below a certain energy 
scale, the so-called Kondo temperature Tk- It is ade- 
quate to ask the question what is the typical length scale 
in which the impurity spin gets screened. As a simple es- 
timate, just by comparing the scales of the competing ki- 
netic and binding energies, one easily gets a length scale, 
the so-called Kondo coherence length, = Hvp/ksTK- 
By substituting the typical Fermi velocity vp in metals, 
and taking a typical Kondo temperature Tk ~ IK, one 
gets £k ~ Ifim, which is comparable with, or even larger 
than the typical dimension of today's mesoscopic devices. 
However, to measure that screening cloud is highly non- 
trivial from experimental point of view. 

In order to observe the screening cloud, one has to mea- 
sure correlations between the impurity spin and the con- 
duction electron spin density. It is very difficult to imag- 
ine such kind of a measurement in bulk metallic samples. 
Recently, it has become possible to perform scanning 
tunneling microscopy (STM) measurements on metallic 
surfaces addressing a single magnetic impurity and its 
close neighborhood^]. In those experiments the screen- 
ing cloud was not observed since the STM tip probes 
the local charge density of states which is affected by the 
Kondo correlations in the close neighborhood of the im- 



purity only. Even in case of having spin-polarized STM 
tips, it is very difficult to measure the spin correlations 
since the Kondo effect manifests itself in continuous, co- 
herent spin flip processes at a time scale ~ 1 /Tk meaning 
that one has to measure at a frequency of tens of giga- 
hertzs. The idea of suppressing the impurity spin fluc- 
tuations by applying a local magnetic field is doomed to 
fail since a magnetic field which is large enough to sup- 
press the spin fluctuations is necessarily large enough to 
destroy Kondo correlations as well. 

There were experimental setups proposed in order to 
measure the Kondo cloud in confined systems^, e.g. in 
case of a quantum dot attached to one dimensional lead 
with length shorter than, or comparable with £k- Such 
kind of a proposal suffers from the facts that it is de- 
structive, i.e. the Kondo effect disappears for one di- 
mensional leads shorter than £jj and one can still argue 
that the suppression of the Kondo effect is a result of 
having the level spacing in the lead 8e > Tk for short 
enough leads therefore one does not need the theoretical 
construction of the screening cloud to explain the exper- 
imental findings. In case of an impurity embedded into a 
higher dimensional environment, another length scale Ik 
emergesQ by equating A(Z#) = Tk where A(Ik) stands 
for the mean level spacing in a box with size Ik- For 
the case of D > 2 that length scale can be substantially 
smaller than fg-. In ID systems, which are in the focus 
of the present paper, the two length scales are essentially 
equal. 

Despite those challenges, there were proposals pub- 
lished to observe the Kondo screening cloud. These pro- 
posals deal with the Knight shift persistent current Q 
or conductance of mesoscopic systems (lpj. 

Very recently, Hand and his coworkers established a 
proportionality between the weight of Kondo resonance 
and the spatial extension of the Kondo correlations in 
mesoscopic systems [H|]. This fact allows another, spec- 
troscopic way to observe the Kondo screening cloud. 
Those proposals are very promising even though the ex- 
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perimental realizations are yet to come. 

In the present paper we are less ambitious since our 
goal is not to propose new ways of observation but to 
extend the powerful and numerically exact NRGfl2| to 
be capable to compute spatial correlations and therefore 
make predictions for the outcome of future experiments. 

Model and definition of correlation functions — The 
Hamiltonian of the Kondo model can be written as 



(i) 



where S is the impurity spin (S = 1/2, located at the ori- 
gin), c\. (c fe „) creates (annihilates) a conduction electron 
with momentum k and spin /i, while <?(0) = ifi (0)c?V>(0) 
stands for the electron spin density at the position of the 
impurity. As it is mentioned above, the system described 
by Eq.JT]) forms a singlet below the Kondo temperature 
Tk ~ Do exp (— 1/2qqJ), where Dq is the high energy 
cutoff and go is the density of states of the conduction 
electrons. (For sake of simplicity, we consider a band 
with a constant density of states in the following. More- 
over, we restrict ourselves for ID electrons because it is 
easier to treat them numerically. When it is adequate, 
we will comment the case of higher dimensions as well.) 
A quantity which measures the spatial extension of the 
singlet is the equal-time correlator of the impurity spin 
and conduction electron spin density at position x: 



Xt=o{x) = (Sa(x)) t =o- 



(2) 



Since we know that the ground state of the Hamiltonian 
Eq.(Q} is a singlet, we can easily derive a sum rule for the 
equal time correlator at zero temperature: 



(Sa(x))t=odx 
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(3) 



Method — To calculate spatial correlations is a non- 
trivial task since most of the methods used to investigate 
the Kondo model are not able to reproduce correlation 
functions. Comparatively only a very little of the the- 
oretical study has been focused on spatial correlations: 
Perturbative calculations have been performed 13J, as 
well as Monte Carlo analysis [15j]- ln the present paper 
we extend NRG to compute the Kondo correlations. 

Wilson's NRG -the most useful numerically exact 
method to obtain correlation functions suffers from 
the very bad spatial resolution away from the impu- 
rity. This fact is a direct consequence of the cor- 
nerstone of the method, the logarithmic discretization 
of the conduction band. In Wilson's NRG technique 
one introduces a discretization parameter, A and with 
the help of that discretizes the conduction band loga- 
rithmically by dividing it into intervals (i.e. the nth 
interval is ] — DqA^"; — DoA~ n ~ 1 ] for negative, and 
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FIG. 1: (a) Spherical shells in real space depicting the extents 
of wave functions of states which are represented as on-site 
states on the Wilson chain. As shown, the NRG method has 
a good spatial resolution around the impurity only, (b) A 
straightforward extension of NRG tackles that problem pro- 
viding a good spatial resolution at the impurity site as well 
as at another, freely chosen point by introducing a second 
electron field. 



[.DoA - ™ -1 ; DoA~™[ for positive energies measured from 
the Fermi energy and Dq stands for the half bandwidth) 
and keeping one mode per interval only. As a next 
step, one maps the problem onto a semi-infinite chain 
with the impurity at the end by means of a Lanczos- 
transformation. As a consequence of the logarithmic dis- 
cretization, the hopping amplitude falls off exponentially 
along the chain allowing us to diagonalize the problem 
iteratively. As it is shown in Ref.[12( the states repre- 
sented as on-site states on the Wilson-chain correspond 
to extended states in real space with a typical spatial 
extension 
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where N is the site index along the chain. It is obvious 
that the numerical renormalization group method has a 
good spatial resolution at the position of the impurity 
only. [See FigHJi.] 

To tackle that problem, we introduce two electron 
fields, one is located at the impurity site, -00 = "0(0) and 
another one at point x where we want to evaluate the cor- 
relation function, ip± = ^>(x).[See FigQJ).] (A very similar 
procedure was applied in the numerical renormalization 
group treatment of the two impurity problem (l7|.) Of 
course, these fields are not orthogonal, i.e. they do not 
fulfill the canonical anticommutational relations. As a 
next step, we can introduce the proper linear combina- 
tions of these fields, 



rl>± = l/V2(Vo±V>i) 



(5) 



which are now anticommuting but the corresponding den- 
sity of states is modified i.e. it acquires energy de- 
pendence, e.g. for ID electrons g]P(s) — Qo/2[l ± 
cos(ex/vf)]- (In 2D the density of states reads e± D (e) = 
Qo/2[l ± Jo(er/vF)] where Jo is the zeroth Bessel func- 
tion, while in 3D Qj^ie) — go/2[l±sin(er/uj?)/(er/ui?)]. 
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FIG. 2: (Color online) The equal time spin-spin correlation 
function, x^=o(x) = (Sa(x)) as a function of the distance 
x measured from the impurity for different values of the 
Kondo coupling j — qqJ. xt=o{x) oscillates as ~ cos 2 (fc_F2;). 
As shown in the inset the envelope of the oscillating part 
(i.e.xt=o(a; = nn/kp) where n is a integer) crosses over from 
~ s _1 to ~ x~ 2 at around the Kondo coherence length, £k- 
The envelope function for different couplings nicely collapse 
into one universal curve. 



As we see, the oscillations of the even/odd density of 
states decay as ~ r ^(- D - 1 )/ 2 j where D is the dimension- 
ality of the problem.) In the new basis the Kondo Hamil- 
tonian reads 



H = 



E 

i=±M=U 



(6) 



In this representation 3{x) = l/2[^ + — ijjl_}<f[ip+ — ip-] 
is represented in a high accuracy. In brief, the main idea 
behind the above transformation is that one can get a 
very good spatial resolution at the impurity site as well 
as at another freely chosen position x if one is willing to 
pay the price of (a) having two electron channels mixed 
by the interaction, (b) having energy dependent density 
of states and (c) having a separate NRG iteration for 
every different x values. Those difficulties are possible to 
be handled: The most serious one amongst them is the 
extension of NRG to arbitrary density of states which has 
been solved by Bulla et al.[18| 

Results — In this NRG scheme the calculation of the 
equal time correlation function is rather simple: Xt=o{ x ) 
appears to be a static thermodynamic quantity which can 
be evaluated with a high precision. The results for differ- 
ent Kondo couplings are shown in Figj2] the correlation 
function oscillates as ~ cos 2 (kpx). As it is shown in the 
inset of Fig|2j the envelope of the correlation function for 
different Kondo couplings nicely collapse into one univer- 
sal curve (apart from the points x < it/hp which show 
non-universal, coupling-dependent behavior, not plotted 




FIG. 3: (Color online) The envelope of Xt=o(x) for different 
couplings and different temperatures: Any finite temperature 
introduces another energy scale, £t = Uvp/kBT at which the 
envelope function crosses over from an algebraic to an expo- 
nential decay. 



in the inset). The envelope curve (i.e.xt=o{% = nit/kp) 
where n is an integer) changes its decay from ~ l/x to 
~ l/x 2 at around Our results supports the heuris- 
tic picture of the Kondo ground state: The impurity is 
locked into a singlet with an electron whose wave func- 
tion is spread out over a distance of If we test the 
correlations on a length scale much shorter than £k the 
impurity appears to be unscreened and a perturbative 
analysis can be performed yielding Xt=o ~ \/x. On 
the other hand, the impurity is screened for distances 
x 3> S,k and Nozieres' Fermi liquid theory applies which 
predicts ~ l/x 2 decay in 1D[19(. There is, however, a 
rather wide crossover regime for which there is no reli- 
able analytic method of calculating Xt=o{%)- From per- 
turbative aspect this is the region where the logarithmic 
corrections to ~ l/x becomes large [H]]. Alternatively, 
following the argumentation of the authors of Refflll. this 
is the region where the impurity can be described by a 
fluctuating spin resulting in exponential decay which is 
overriden by the ~ l/x 2 decay for even larger distances. 

To plot the universal curve in the inset of Fig[2] we 
have chosen the definition of Tk based on the NRG flow 
of the energy levels: We identified Tk as the energy scale 
at which the energy of the first excited state reaches 80% 
of its fixed point value. We have taken this definition 
since the energy spectrum is the most basic and most 
reliable result of any NRG calculation and -as a conse- 
quence of the universality of the Kondo effect- the value 
of Tk obtained in this way can differ from that of other 
definitions up to a constant multiplying factor only. 

Note that the results shown in Figj2] correspond to ID: 
in the case of higher dimensions the correlation function 
tends to zero even faster since the sum rule (Eq[3]) fixes 
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its integral. This fact enters the numerical calculation 
through the ~ r ^(- D_1 )/ 2 decay of density of states oscil- 
lations and makes the calculation of Xt=o( x ) challenging 
for dimensions higher than D = 1. 

Up to this point we have considered the case of zero 
temperature T = 0. However, it is known that the tem- 
perature plays an essential role in Kondo physics: Kondo 
correlations are destroyed when the temperature reaches 
the Kondo scale T ~ Tr- It is adequate to ask the ques- 
tion how does this fact show up in Xt=o(x)t 

Since Xt=o(x) is a static quantity we obtain its com- 
plete temperature dependence by numerical renormaliza- 
tion as the iteration proceeds. In Fig|3] we show the re- 
sults for T = 0, T = T K /Z, T = T K , and T = 3T K . As it 
is transparent from the curves, the effect of finite temper- 
ature shows up in the appearance of another length scale, 
the thermal length scale, £t = hvF /ksT. For x < £t the 
correlation function Xt=o( x ) is not much affected while 
for x > I^t the correlations are cut off exponentially 

Based on those results we can now interpret the role 
of finite temperature in Kondo screening as follows: The 
fact that the impurity spin is perfectly screened at T = 
is reflected in the sum rule given by Eq.Q. At any fi- 
nite temperature the integral is reduced and can approx- 
imately rewritten as J^ T Xt=o(x)dx. When the temper- 
ature is lower than the Kondo temperature T < Tk-, 
the corresponding thermal length scale is larger than £k 
Since Xt=a( x ) nas its most weight in the region x < 
the effect of a small temperature is just a small correc- 
tion to the perfect screening obtained at T = 0. Given 
that in the x > £k regime Xt=o(x) ~ x~ 2 , the correction 
appears to be linear in T. In contrast, for T > Tk the 
corresponding £t is shorter than and the correction 
to the integral is not small and consequently, we cannot 
speak about the screening of the local moment any more. 

Conclusions — We have shown that Wilson's NRG 
technique is capable to handle spatial correlations if we 
apply a straightforward extension. To demonstrate that 
we have computed the (Sa(x))t=o correlator for a one 
dimensional Kondo system and shown that the decay of 
spin correlations crosses over from ~ x _1 to ~ x~ 2 at 
around the Kondo coherence length, = hvp/ksTK- 
We have shown by calculating the temperature depen- 
dence that any finite temperature introduces a new en- 
ergy scale beyond which the Kondo correlations vanish 
exponentially Our method is -apart from the numerical 
challenges- easy to generalize to other impurity mod- 
els. A very attractive example would be the two channel 
Kondo model to see how the non-Fermi liquid nature of 
the ground state shows up in the spatial spin correlations. 
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